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Abstract 

The ANTARES detector was operated in a configuration with 5 lines for a period of 10 months from February until 
November 2007. The duty cycle was better than 80% during this period and almost 2 ■ 10^ atmospheric muon triggers 
were collected. This large sample was used to test Monte Carlo simulation programs and to evaluate possible systematic 
effects due to uncertainties on environmental parameters and detector description. First results are presented and 
discussed. 
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1. Introduction 

The ANTARES Collaboration has just com- 
pleted the construction of a neutrino telescope off 
the Southern French coast, 40 km from Toulon, at 
a depth of about 2500 m under the sea level. Tech- 
nical details on the ANTARES detector and first 
results are presented and discussed at this workshop 
[T]. The main goal of ANTARES is the detection 
of high energy neutrinos of astrophysical origin. 
Nevertheless, the most abundant signal is due to 
high energy muons remaining from the extensive 
air showers produced in interactions between pri- 
mary cosmic radiation and atmospheric nuclei. The 
shielding effect of the sea reduces their flux, which, 
however, is several orders of magnitude larger than 
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the atmospheric neutrino flux. They represent a 
dangerous background for track reconstruction as 
their Cherenkov light can mimic fake upward going 
tracks. On the other hand they are a useful tool to 
test offline analysis software, to check our under- 
standing of the detector and to estimate systematic 
uncertainties. 

From February to November 2007, the detector ran 
in a 5- line conflguration and collected almost 2 • 10^ 
events mainly due to atmospheric muons. This large 
amount of data was used for comparisons with the 
results from a Monte Carlo (MC) simulation. More- 
over, several MC samples were produced, using dif- 
ferent sets of input parameters, to give an estimate 
of the systematic effects due to environmental and 
geometrical parameter uncertainties. 



2. Monte Carlo simulation 

Two different approaches can be adopted to sim- 
ulate the atmospheric muon flux: 
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- a parameterized description of underwater muon 
flux, 

- a full Monte Carlo simulation. 

In ANTARES software both methods are used. 
The former is based on the MUPAGE program, de- 
scribed at this conference in |:2 . 
In the present analysis the second approach is used 
for comparison to data and to evaluate systematic 
errors. In the following its general scheme is pre- 
sented. 

The first step of the full MC simulation consists of 
the generation of a large number of air showers in- 
duced by primary nuclei with energy ranging from 
1 to 10^ TeV/nucleon and zenith angles between 0° 
and 85" using the CORSIKA software (vrs.6.2) [3] 
and the hadronic interaction model QGSJET.Olc 
The energy spectrum used for generation is 
E~'^. Then, all muons reaching the sea level, with 
energy larger than a threshold energy, are prop- 
agated through sea water to the detector, using 
MUSIC (MUon Simulation Code) a 3D muon 
propagation program that takes into account the 
main energy loss mechanisms for muons. 
Finally, muons are transported through the 
ANTARES sensitive volume, Cherenkov light is 
produced and the detector response is simulated. 
At this point, files containing information on hit 
photomultipliers (PMTs) are translated into the 
same format used for data, after the addition of a 
background chosen by the user and a trigger se- 
lection, and are ready to be treated with the same 
reconstruction and analysis program chain used for 
data. The composition of the primary cosmic ray 
fiux is defined by the user. In this study, a simplified 
version (only 5 mass groups are considered) of the 
model presented in [6j (Horandel model) is used. 

3. Systematic effects 

During MC simulation several input parameters 
are required to define the environmental and ge- 
ometrical characteristics of the detector. Some of 
them play a role as sources of systematic uncertain- 
ties. 

In the following the effect of water absorption length 
and of PMT description on the reconstructed track 
rate is considered. 

The absorption length of light in water depends on 
the light wavelength and has been measured during 
several sea campaigns. The reference values used in 
the MC simulation are the result of a fit to a set of 
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Fig. 1. Absorption length used as input for ANTARES MC 
simulation. The lower line is a fit to measurements (black 
points) taken in the ANTARES site. The upper line is the 
absorption length for pure water and is shown for compari- 
son. 

measurements collected in the ANTARES site, fig. 
[l] Varying by ±10% the reference values of absorp- 
tion length, 2 new samples of atmospheric muons 
were created and compared to the standard produc- 
tion. The 3 zenith distributions are shown in fig. [2j 
There is an almost negligible effect on the shape of 
the histograms, while the absolute fiux changes by 
+25%/ -20%. 

The geometrical description of the photomultiplier 
implies the knowledge of its effective area, which 
is the result of a convolution of several factors like 
quantum efficiency, geometrical surface of the pho- 
tocathode and collection efficiency. The standard 
MC simulation takes as reference values what is re- 
ported on Hamamatsu technical sheets. A sample of 
simulated muons has been prepared using a PMT 
efficiency decreased by 10%, considering the official 
Hamamatsu values as the most optimistic set of pa- 
rameters. A decrease of about 15% was observed in 
the muon ffux. Decreasing the effective area by the 
same amount produces a reduction of 7% of the at- 
mospheric neutrino ffux. 

Finally, the effect of the maximum angle between the 
PMT axis and the Cherenkov photon direction (t?c, 
see fig.|3| allowing light collection has been studied. 
In fig. ^the OM efliciency is shown as a function 
of dc- Presently, a cut-off at cos(i9c)=-0.8 is used. It 
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Fig. 2. Zenith distributions of reconstructed tracks using dif- 
ferent absorption lengths. The upper (lower) histogram cor- 
responds to an increase (decrease) of the absorption length 
by 10 % wrt the reference values. 




Fig. 3. Description of the angle i?c used to define the limit of 
detection of Chorenkov photons. Two cases corresponding to 
a downward going and to an upward going muon are shown. 

means that for -de larger than ^ 140° no photon de- 
tection is expected. Moving the cut-ofFto -0.65 (cor- 
responding to iDc ~ 130°) a decrease of the rate of 
reconstructed tracks of about 20% is observed, see 

fig-Hi 

Given PMTs' orientation (looking downward at 45° 
with respect to the vertical) the effect on neutrino 
detection is lower. 

Summing in quadrature the different contributions, 
a global systematic effect of about ±30% can be 
considered as an estimate of the errors produced by 
uncertainties on environmental and geometrical pa- 
rameters. New measurements are in progress both 
of absorption length, using in-situ tools, and of ge- 




Fig. 4. Parameterization of OM angular efficiency with the 
indication of the two different cutoff values used in the sys- 
tematic effect study. 
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Fig. 5. Effect of the angle i?c cut-off on zenith distribu- 
tion of reconstructed tracks: solid line refers to cos()9c)=-0.8 
(reference value in MC simulation), dotted line represents 
cos (i?c) =-0.65. 

ometrical characteristics of the PMTs and a more 
precise definition of the actual values is expected. 
Another source of uncertainty, which is common to 
all astroparticle physics experiments, is represented 
by the choice of the primary cosmic ray composition 
and of hadronic interaction description. Some ex- 
perimental measurements suggest that the Horandel 
model [6], adopted in this MC simulation, underes- 
timates the all-particle spectrum by about 30% (see, 
for example [7] and [5]). 

A similar uncertainty can be considered as an esti- 
mate of the systematic effect due to the choice of the 
hadronic interaction model. 
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Fig. 6. Zenith distribution of reconstructed tracks. Black 
points represent data. Lines refer to MC expectations: solid 
line = full simulation with CORSIKA code + QGSJETOl 
for hadronic interaction description and Horandel model, 
[6]; dotted line = full simulation with CORSIKA code + 
QGSJETOl for hadronic interaction description and NSU 
model, |9]; dashed-dotted line = MUPAGE parameteriza- 
tion, [2]. The shadowed band represents the systematic er- 
ror, starting from the Horandel model, due to environmental 
and geometrical parameters. See text for details. 

To give an idea of the possible difference among 
different parameterizations of primary cosmic ray 
composition, in fig.|6]the zenith distribution of data 
(black points) is compared to MC expectations ob- 
tained using the Horandel model (solid line) [6J , the 
NSU model (dotted line) 9^ and the MUPAGE pa- 
rameterization (dashed-dotted line) 2 . The shad- 
owed band represents the systematic error due to 
environmental and geometrical parameters. In fig. 
[Tjthe azimuth distribution is shown. Only data and 
MC with Horandel model are shown, together with 
the systematic error band. The pictures show a good 
agreement between data and MC both in shape and 
in absolute normalization, within the estimated un- 
certainties. 

4. Conclusions 

A comparison between data and results from MC 
simulation of atmospheric muon flux is presented. 
The distributions of zenith and azimuth show a 
good agreement both in shape and in the absolute 
flux. 

The effect of environmental and geometrical pa- 
rameter uncertainties on reconstructed track rate 
is studied. Presently, a value of about ±30% can be 
considered as a first evaluation of the systematic 
errors. New measurements are in progress to reduce 
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Fig. 7. Azimuth distribution of reconstructed tracks. Black 
points represent data. The line refers to MC expectations 
obtained from the full simulation with CORSIKA code -|- 
QGSJETOl for hadronic interaction description and Horan- 
del model. The shadowed band represents the systematic 
error due to environmental and geometrical parameters. See 
text for details. 

uncertainties on the parameters used as input for 
MC simulations. The systematic effect on atmo- 
spheric neutrino flux predictions is noticeably lower. 
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